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' Abstract 

o 

. The large-/?o limit of QCD is discussed, with the emphasize on simple technical 

methods of calculating various quantities at the order 1/Po- Many examples, mainly 
from heavy quark physics, are considered. Some QCD results based on renormaliza- 
C*~) ■ tion group (and not restricted to the large-/?o limit) are also discussed. 

O 

1 Large-/3o Limit 

It is well known that perturbative series do not converge. They are asymptotic series, i.e., 
the difference between the exact result and its approximation up to the order af, divided 



by , tends to in the limit a s — ► 0. Large-order behaviour of various perturbative series 
attracted considerable attention during recent years. Most of the results obtained so far 
are model-dependent: they are derived in the large-/3o limit, i.e., at rif — > — oo. There are 
some hints that the situation in the real QCD may be not too different from this limit, but 
this cannot be proved. However, a few results are rigorous consequences of QCD; they are 
based on the renormalization group. Many more applications are discussed in the excellent 
review jTj, where additional references can be found. 

Let's consider a perturbative quantity A such that the tree diagram for it contains 
no gluon propagators. We can always normalize the tree value of A to be 1. Then the 
perturbative series for the bare quantity Aq has the form 



oo L— 1 / 

( 4vr )d/5 



L=l n=0 



^ = l + EE^/(77^) , (1-1) 
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where L is the number of loops. This series can be rewritten in terms of /3 = yC^ — f 
instead of nf 




oo L— 1 




) 



L 



(1.2) 



L=l n=0 



Now we are going to consider f3 as a large parameter such that (3 a s ~ 1, and consider 
only a few terms in the expansion in l//5 ~ ot s : 



This regime is called large-/?o limit, and can only hold in QCD with nf — > — oo. Note that 
it has nothing in common with the large- N c limit, because we cannot control powers of N c 
in the coefficients am- 

There is some empirical evidence [2] that the two-loop coefficients a2i/?o + 0,20 for many 
quantities are well approximated by a2i/?o- It is easy to find a 21 from the diagram with the 
quark-loop insertion into the gluon propagator in the one-loop correction. Then we can 
estimate the full two-loop coefficient as a 2 i/3o — a 2i { n f ~~ ^Ca/T f ). This is called naive 
nonabelianization (2]. Of course, there is no guarantee that this will hold at higher orders. 
We can only hope that higher perturbative corrections are mainly due to the running of 
a s ; in this respect, gluonic contributions behave as —33/2 flavours, and QCD with nf = 3 
or 4 flavours is not too different from QCD with —00 flavours. 

It is easy to find the coefficients oll,l-i of the highest degree (3q _1 at L loops. They are 
determined by the coefficients a' L L _ x of i-e., by inserting L — 1 quark loops into the 

gluon propagator in the one-loop correction. We shall assume for now that there is only 
one gluon propagator and no three-gluon vertices at one loop. The bare gluon propagator 
with L — 1 quark loops inserted is 



in the Landau gauge ao = 0, with a shifted power of — p 2 (and an extra constant factor). 
Let the one-loop contribution to A be (a x + a[a ) g 2 /(An) d l 2 . If we calculate the one-loop 
contribution in the Landau gauge a\ with the denominator of the gluon propagator equal 
to (— p 2 ) n instead of just — p 2 and call it ai(n), then for all L > 1 






(1.4) 



It looks like the free propagator 
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Only the one-loop contribution a w contains the additional gauge- dependent term a[a . 

The large-/5o limit, as formulated above, does not correspond to summation of any 
subset of diagrams. If we include not only quark loops, but also gluon and ghost ones, 
then —-sTprif in (jl.4|) is replaced by 



C A L 3-2e 



3 [ 2(1 



Po-^r{8e + — r(a + 3) l--(a + 3) }. (1.6) 
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Summing these diagrams yields a gauge- dependent result. This gauge dependence is com- 
pensated (for a gauge-invariant A ) by other diagrams, which have more complicated 
topologies than a simple chain, and are impossible to sum. In the gauge a Q = —3, one-loop 
running of a s is produced by one-loop insertions in the gluon propagator only, without ver- 
tex contributions. Summation of chains of one-loop insertions into the gluon propagator 
in this gauge is equivalent to the large-/3 limit. 

In the large-/9 limit, Pi ~ p , P2 ~ Po, etc. Therefore, /5-function is equal to 

0=^ (1.7) 

(this term is of order 1) plus 0(1/ P ) corrections. At the leading order in 1/ /? , the 
renormalization-group equation 



d log Z a (3 



dp e + p 

can be explicitly integrated: 
1 



1.9) 



" l + P/e' 
To the leading order in l//3o, 

-M = Afe- (L10) 

'His 

The perturbative series ()1.2|) can be rewritten (in the Landau gauge) via the renormal- 
ized quantities: 



00 



+ OI — I, (in! 



L=l 

where 

F{e, u) = ue 7£ ai(l +u- e)fx 2u D(s) u/£ - 1 . (1.12) 

If a 7^ 0, the term a[aP/Po + (9(1//3q) should be added (the difference between ao and a is 
O{l/Po)). 
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We can expand in the renormalized a s , or in (3 fll.7|) . using 



P V f^ L 



e + (3 



(here (x) n = x(x + 1) • - • (x + n — 1) = T(x + n) /T(x) is the Pochhammer symbol). In the 
applications we shall consider, F(e, u) is regular at the origin: 

oo oo 

F(e,u) = Y,Ys F ™ £num > ( L13 ) 

n=0 m=0 

though I know no general proof of this fact. Substituting these expansions to (jl.llj) . we 
obtain a quadruple sum expressing A via the renormalized quantities. 

The bare quantity A = ZA, where both Z and A have the form l + (9(l//?o). Therefore, 
we can find Z — 1 with the l/flo accuracy just by retaining all terms with negative powers 
of £ in this quadruple sum. The renormalized A — 1, with the l//?o accuracy, is given by 
terms with e°. It is enough to find Z\, the coefficient of 1/e in Z, in order to have the 
anomalous dimension 

dZ l 

7 = -2- 



d\og(3 

Collecting terms with e~ l in the quadruple sum for A , we obtain for /3qZi 

(3F 00 - (3 2 {F 10 + F i) + (3 3 {F 20 + F u + F 02 ) - /? 4 (F 30 + F 21 + F 12 + F 03 ) + • • • 
+ \P 2 {F XQ + 2F i) - (3 3 (F 20 + 2F n + 4F 02 ) + §/3 4 (F 30 + 2F 21 + 4F 12 + 8F 03 ) + 
+ l(3 3 (F 20 + 3F U + 9F 02 ) - /3 4 (F 30 + 3F 21 + 9F 12 + 27F 03 ) + • ■ ■ 
+ \f3\Fw + 4F 21 + 16F 12 + 64F 03 ) + • • • 
_l_ . . . 



(3 2 p 3 (3 A 
= fiFoo — 2"^ 10 + ~3"^ 20 — 4~ F30 ~^ 

Therefore, the anomalous dimension is |3j 

P , „ ( 1 



Po K J \p t 



7 = -2^F(-/3,0) + O . (1.14) 



Collecting terms with e° in the quadruple sum for A Q , we obtain for Pq(A — 1) 

P(F 10 + Foi) - P 2 (F 20 + F n + F 02 ) + /? 3 (F 30 + F 21 + F 12 + F 03 ) 

— /3 4 (F 40 + F 3 i + F 22 + F13 + F04) + ■ ■ • 
+ ^ 2 (F 20 + 2F U + 4F 02 ) - p 3 (F 30 + 2F 21 + 4F 12 + 8F 03 ) 
+ |/? 4 (F 40 + 2F 31 + 4F 22 + 8F 13 + 16F 04 ) + • • • 
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+ §/3 3 (F 30 + 3F 21 + 9F 12 + 27F 03 ) - (3\F m + 3F 31 + 9F 22 + 27F 13 + 81F 04 ) + • • • 
+ \f3\Fw + 4F 31 + 16F 22 + 64F 13 + 256F 04 ) + • - - 
+ ... 

(3 2 (3 3 (5 4 
= (3F 10 — ^"-^20 + "3" "^30 — 4" "^40 H 

+ /3F 01 + /3 2 F 02 + 2/? 3 F 03 + 6/? 4 F 04 + • • • 

Therefore, the renormalized quantity is jl] 

00 

m = i + i l /(m)-f(o.q) 1 f fe ^M)-fW) +0 m _ (L15) 

/W ^ pW « \/3o/ 

-/3 

where (3 = j3oa s (/j,) / (iTr) . 

The renormalization group equation 

d log A(fi) j(a s ) 



dloga s 2(3(a s ) 
can be conveniently solved as 

A{y) = A ( -^4 ) ^ K 7 (a s (/i)) , (1.16) 



where the function 

~ , , f ( i{ol s ) 70 \da s 70 / 71 (3A a s 

K lM = ex P y ^ - -j - = 1 + ^ [- - ^ + - (1.17) 



satisfying 

-ftTo(as) = 1 , if_ 7 (a 5 ) = iif~ 1 (a a ) , K J1+J2 (a s ) = K 7l (a s )K l2 (a s ) 
has been introduced, and 

A = A(fi )K^(a s (fx )) . 
At the first order in 1/po, we obtain from (|1.14jl 

1 f . F{e, 0) - F(0, 0) 



if 7 (a a ) = 1 + — / cfe 
Po J 

Therefore, 



00 

1 r^„-uM« s u 0) ) F(0,u)-F(0,0) ,1 1 



i = 1 + — / due- ulp{ - aa ^ o)) 



Po J u 



+ 



ft 



Let's suppose that m ^> A^jg is the characteristic hard scale in the quantity A. Then 
F(e,u) contains the factor (fx/m) 2u . When taking the limit e — > 0, the factor D{e) u ^ e ~ l 
in 1)1.12)1 becomes exp (§«)■ Therefore, 

F(Q,u)=(?—^] F(u), F(u) =ua 1 (l + u)m 2u . (1.18) 



m 

It is most convenient to use 

^ = e -B/fl m (1.19) 

in the definition of A. In the rest of this Chapter, (3 will mean f3 a s (no) / (47i) . This 
renormalizat ion-group invariant is 

00 

A = 1 + i- / d M e- u//3 5(M) + J > ( L2 °) 



where 



F M -F(0) 

It 



Here, 



; - if ' 5/6 ^MS 



m 



If we substitute the expansion 

00 

S(u) = ]T s L u L ~ x 

L=l 

into the Laplace integral (|1.2()j) . we obtain the renormalized perturbative series 

i. 1 + _L| c ^ + (_L), (1 . 23) 



1.24) 



^ = ( L - 1)!s -(s)"' SM 

u=0 

Therefore, S(u) can be obtained from A (jl.2Hj) by 

00 

which is called Borel transform. 

We see that the function F(e, u) (Jl. 12)1 contains all the necessary information about the 
quantity A at the 1/po order. The anomalous dimension 1)1.14)1 is determined by F(e,0), 
and the renormalization-group invariant A (Jl. 20)) (which gives (Jl.lfijl ) - by F(0,u). 

These formulae are written in the Landau gauge a = 0; if a 7^ 0, additional one- loop terms 
from the longitudinal part of the gluon propagator should be added. 
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2 Renormalons 



The Laplace integral (jl.2()j) is not well-defined if the Borel image S(u) has singularities 
on the integration path - the positive half-axis u > 0. At the first order in 1//3q, S(u) 
typically has simple poles. If 

S(u) = — !— + ... (2.1) 

Uo — u 

where dots mean terms regular at u = u , and u > 0, then the integral (jl.20j) is not 
well-defined near uq. One way to make sense of this integral is to use its principal value: 
to make a hole [uq — 5, uq + S] and take the limit S — > 0. However, if we make, e.g., a hole 
[uq — 5, uo + 25} instead, we'll get a result which differs from the principal value by the 
residue of the integrand times log 2. Therefore, the sum of the perturbative series ()1.20j) 
contains an intrinsic ambiguity of the order of this residue. It is equal to 

rp -u //3 / P 5/6A \ 2U ° 

^= r ^— = ^( e —^) . (2.2) 
Po Po V m J 

These renormalon ambiguities are commensurate with 1/m power corrections - contribu- 
tions of matrix elements of higher-dimensional operators to the quantity A. The full result 
for the physical quantity A must be unambiguous. Therefore, if one changes prescription 
for handling the integral across the renormalon singularity at u = Uq, one has to change 
the values of the dimension-2wo matrix elements accordingly. This shows that renormalons 
can only happen at integer and half-integer values of u, corresponding to dimensionalities 
of allowed power corrections. The largest ambiguity is associated with the renormalon 
closest to the origin. 

The renormalon pole (j2.1j) yields the contribution to the coefficients cl of the renor- 
malized perturbative series (jl.2Hj) equal to 

cl = r y —^ 2.3 

(see (jl.24j) ). The series (j2.1|) is, clearly, divergent. Using the Stirling formula for the 
factorial, we can see that the terms of this series behave as 



eu Q 

at large L. The best one can do with such a series is to sum it until its minimum term, 
and to assign it an ambiguity of the order of this minimum term. The minimum happens 
at L fa Uq/(3 loops, and the magnitude of the minimum term is given by (|2.2j) . This is 
another way to look at this renormalon ambiguity. The fastest-growing contribution to cl 
comes from the renormalon most close to the origin. 

Note that renormalons at Mo < give sign-alternating factorially-growing coefficients (|2.Hj) . 
For such series, the integral (jl.2()|) provides an unambiguous definition of summation called 
Borel sum. 
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Renormalon singularities can result from either UV or IR divergences of the one-loop 
integral. Suppose that it behaves at k — > oo as j d A k/(—k 2 ) nvv , so that the degree of 
its UV divergence (at d — 4) is i^uv = 4 — 2nuv- When we insert the renormalon chain, 
the power changes: riuv — > ^uv + (L — l)e = ^uv + u if £ = (which is the case when 
calculating S(u)). This integral can have an UV divergence only at u < 2 — ri\jy = u\jy/2. 
Therefore, UV renormalons can be situated at z^uv/2 and to the left. Only quantities 
A with power-like UV divergences at one loop have UV renormalons at positive u. The 
divergence at u = is the usual UV divergence of the one-loop integral, which is eliminated 
by renormalization; renormalized quantities have no UV renormalon at u = 0. 

Similarly, if the one-loop integral behaves as J d A k/(—k 2 ) nm at k — ► (where k is the 
virtual gluon momentum), so that the degree of its IR divergence is Vm = 2niR — 4, S(u) 
can have an IR divergence only at u > 2 — n = —vir/2, and IR renormalons can be situated 
at — vm/2 and integer and half-integer points to the right from it. Quantities described by 
off-shell diagrams have niR = 1, and their IR renormalons are at u — 1 and to the right. 

We can get a better understanding of the physical meaning of renormalons if we 
rewrite (jl.2U|) in the form |3j 







d Z w{T) <^ +0 ri 

r 4tt V/?o 



A = l+ l -w{r) ^ A _^ +0[±- 2 ) . (2.4) 



This looks like the one-loop correction, but with the running a s under the integral sign. 
The function w(t) has the meaning of the distribution function in gluon virtualities in the 
one-loop diagram; it is normalized to the coefficient of a s /(4ir) in the one-loop correction. 
Inside the l//3 term in (j2.4j) . we may use the leading-order formula for running of a s : 

«.(V^A*o) = 1 I'jff = a.wf;(-/31ogrr. 
1 + p log r n=Q 

Substituting this expansion into (|2.4jl . we see, that this representation holds if w(r) is 
related to the perturbative series coefficients cl by 

oo 

c l = I yw(r)(-logr) L - 1 . (2.5) 
o 

Therefore, S(u) (|1.25j) becomes 

oo 

f dr 

S(u) = / — w{t)t' u . (2.6) 
o 

In other words, S(u) is the Mellin transform of w(t). Therefore, the distribution function 
w(t) is given by the inverse Mellin transform: 

uo+ioo 

w(t) = ^- [ duS(u)r\ (2.7) 
2m J 

UQ—ioo 
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where Uq should lie in the gap between IR and UV renormalons. 

At r < 1 we can close the integration contour to the right. If S(u) has IR renormalons 
fi/ (ui — u), then 

W ( T ) = ^V;t Ui . (2.8) 

IR 

The leading term at small r is given by the leftmost IR renormalon. If our quantity A is 
IR-finite at one loop, all Ui > 0, and w(t) — > at r — > 0. Similarly, at r > 1 we can close 
the contour to the left. If the UV renormalons are Tij \u — ui), then 

w(t) = ^r;r" ! . (2.9) 
uv 

The leading term at large r is given by the rightmost UV renormalon. If A is UV-finite at 
one loop, all u^ < 0, and w(r) — > at r — > oo. 

All virtualities (including small ones) contribute to (|2.4j) . Behaviour of the distribution 
function w(t) in the small- virtuality region r — > is determined by the IR renormalon most 
close to the origin. The integral (|2.4jl is ill-defined, just like the original integral (|1.2()jl . 
The one-loop a s (jl.lOj) becomes infinite at r = (e 5 ^ 6 A^/m) 2 (Landau pole), and we 
integrate across this pole. This happens at small r; substituting the asymptotics ()2.8|) of 
the distribution function at small virtualities, we see that the residue at this pole, given 
by the IR renormalon nearest to the origin, is again equal to f!2.2[) . 



3 Light Quarks 

First, we shall discuss the massless quark propagator at the order 1//3q. The one-loop 
expression for the quark self-energy S(p) in the Landau gauge with the gluon denominator 
raised to the power n = 1 + (L — l)e is 

C F f d d k iTrWdk+firf ( , k^K 



-p 2 J vr<V 2 [—{k + p) 2 } {-k 2 ) n V - k2 
Using the one-loop integrals 

d d k 



iir d/2 {-p 2 ) d/2 - ni - n2 G{n h n2) 



-(k + p) 2 }m[-k 2 }^ 

, . _ Tj-d/2 + ni + n 2 )T(d/2 - m)r(d/2 - n 2 ) 1 ' ' 

we can easily find the function F(e,u) (jl.!2|) . Such functions for all off-shell massless 
quantities have the same T-function structure resulting from (|3.1|) with n 2 = 1 + u — e: 

F(r v) - ( ^ ^1 " + M)r(1 ~ M)r(2 ~ £) DlrY'^N(r v) ft 2^ 

F(£ ' M) "V^J 6 T(2 + u-e)T(3-u-e) D{e) * (e ' tt) - (3 ' 2) 
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The first T-f unction in the numerator, with the positive sign in front of u, comes from the 
first T-function in the numerator of (j3.1j) . with the negative sign in front of d, and its poles 
are UV divergences. The second T-function in the numerator, with the negative sign in 
front of u, comes from the second T-function in the numerator of (|3.1J) . with the positive 
sign in front of d, and its poles are IR divergences. For we obtain 



N(e,u) = -C F (3- 2e) (u - e) 



(3.3) 



At one loop (L — u/e — 1), the Landau-gauge self-energy vanishes; at L — 2, the /?o-term 
in the two-loop result is reproduced. 

The massless-quark propagator S(p) with the accuracy in the Landau gauge is 
equal to 1/f times (fl.llj) . where F(e, u) is given by (|3.2(h ()3.3|) . Terms with negative powers 
of e in its expression via renormalized quantities form the quark-field renormalization 
constant Z q . The anomalous dimension is given by p. 14)1 : 

1= .± Etm +0 (±. 

In the general covariant gauge, the one-loop term proportional to a should be added: 



7g 



2a 



P (1 + 1/3) 



B(2 + (3,2 + (3)T(3 + (3)T(l- (3) 



(3.4) 



^[2a + 3/3(l + |/3-l/3 2 + ...)] 



47T 



This perturbative series for ~f q has the radius of convergence equal to the distance from the 
origin to the nearest singularity, which is situated at (3 — —5/2; in other words, it converges 
at \P\ < 5/2. It reproduces the leading-/? terms in the 2-, 3-, 4-loop results |E1IZ|- 

The renormalized expression for ]t>S(p) is given by f)1.15|) . If we factor out its /x- 
dependence as in (11.16(1 . then the corresponding renormalization-group invariant is given 
by (dHH) with 



S(u) 



1 



N(0,u) 



N(0,0) 



(1 + u)(l -u){2-u) 



3C F 



(l+u)(l-u)(2-u) 



(3.5) 



(here \/—p 2 plays the role of m). The pole at u = — 1 comes from the first T- function 
in the numerator of (|3.2(l . and is an UV renormalon; those at u = 1, 2 come from the 
second T-function, and are IR renormalons (Fig. ^i). We can also see this from the power 
counting (Sect. |2J). The light-quark self-energy seems to have a linear UV divergence. 
However, the leading term of the integrand at k — >• oo, j^/(k 2 ) 2 , yields after integration, 
due to the Lorentz invariance. The actual UV divergence is logarithmic: z/yv = 0, and 
UV renormalons can only be at « < 0. The UV divergence at u = is removed by 
renormalization, and UV renormalons are at u < 0. The index of the IR divergence of the 
self-energy, like that of any off-shell quantity, is i^r = —2, and IR renormalons are at u > 1. 
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Power corrections to the light-quark propagator form an expansion in l/(— p 2 ), therefore, IR 
renormalons can only appear at positive integer values of u. For gauge-invariant quantities, 
the first power correction contains the gluon condensate <G 2 > of dimension 4, and the 
first IR renormalon is at u = 2. The quark propagator is not gauge-invariant, and the 
renormalon at u — 1 is allowed. The virtuality distribution function ()2.7)1 is 

w(r) = -3C F x h _ t 3 ' < 

(Fig.Cb). 




Figure 1: UV renormalons (black squares) and IR renormalons (black circles) in the light- 
quark self-energy (a); the virtuality distribution function (b) 

Now we shall discuss light-quark currents 

jM = Zr r l^)q' Tq , r = 7 ^--- 7 ^. 

Let's calculate their vertex function T(p, 0) up to one loop (Fig. EJ). It is convenient to 
rewrite T as 

r+ + r_. i = — ( r± (A p 



2 V" ~ P 2 
where pT^ = aT^p 1 , a — ±1. Then 

7^7^ = 2ah(d)T a , (3.6) 

where for n antisymmetrized 7-matrices 

h(d)= V (n-i) , V = (-l) n+1 a. (3.7) 



The vertex function T(p, 0) for a Dirac matrix r± is T± ■ T(p 2 ), where the scalar function 
T(p 2 ) can be calculated via h(d), once and for all Dirac matrices (see j2]). 
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a 



b 



Figure 2: Proper vertex T{p,p') of a bilinear quark QCD current 



Calculating the vertex function T(p 2 ) in the Landau gauge with the denominator of the 
gluon propagator raised to the power n = 1 + (L — l)e, we obtain Ql.lljl . ()3.2|) with 



N{e, u) = -C F [2-u-e + 2h(u - h)} . 



(3.1 



For the longitudinal vector current (h = 1 — d/2), the result can be obtained from by 
differentiating S(p 2 ) (Ward identity). The anomalous dimension of the current is 



d log Z Tn 



+ 7<? 



dlogfi 

where the derivative of Zm is given by (jl.l4j) . We arrive at |2| 



-C, 



a. 



(n- 1) (3-n + 2/5) 



3 * 4vr B(2 + (3,2 + p)T(3 + (3)T(1 - (3) 

13n - 35 



-2C F ^(n - 1) 



n — 15 „ 
o 



12 



(3 2 + 



(3.9) 



Naturally, it vanishes at n = 1. This perturbative series converges at \(3\ < 5/2. It 
reproduces the leading powers of (3q in the two- and three-loop results IE] In particular, 
the mass anomalous dimension is 



In 



—Jjo = 2C F 



1 + 1/5 



4tt 5(2 + /?, 2 + /3)r(3 + (3)T(l - (3) 



+ o 



1 



(3.10) 



(it is known at four loops |HJ HD] ) • 

There is a general belief that one may use a naively anticommutating 7^° in open quark 
lines without encountering contradictions, see [TT] . The pseudoscalar currents with 7^°, 
Jac(aO = ^p,Ac(/ i )9o7^ C <L and with the Hooft-Veltman 7^, j H v(/u) = Zp^ v (fi)q' 'yf Y q 
are related to each other by a finite renormalization 



Jac(aO = Zp{ot s {y))jm{v) , ^p(«s) 
Similarly, the axial currents are related by 

3aM = Z A(u s (fij)]£ Y (fi) . 



ZP1 4^ + ZP2 V4vr 



a. 



+ 



(3.11) 



(3.12) 
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The finite renormalization constants Zp^{pc s ) can be obtained from the currents' anoma- 
lous dimensions. Multiplying the Dirac matrix Y in the current by 7^° does not change 
the anomalous dimension; multiplying by 7^ v means n — > 4 — n. Differentiating (|3.11J) 
and ()3.12|) . we have 



d\ogZ P (a s ) _ j j0 (a s ) - j j4 (a s ) 



d log a s 
dlogZ A (a s 
d log a s 

Therefore, 

Z P (a s ) = K. 



(3.13) 



where 7ji = 



"/jo-"/jA a s) J 



Za[ol s 



K, 



(3.14) 



(see (I1.17j) ). For the tensor current, multiplying a^ v by 75 is merely a space-time trans- 

-ia 23 , and thus it does not change the anomalous dimension. 



formation, e.g., 7^ v o~ 01 
Therefore, the constant relating the currents q'j^a^q and q'^ v a^ v q is 

Z T (a s ) = 1 . 

In the large-/3 limit, we obtain from ([3. 13)1 . ()3.9|) 



(3.15) 



1 _4-Cf 

3(3 



1-4CV— 



dp 



Air 



B(2 + P,2 + P)T(3 + p)r(l-p) 



1 + L p _ + 

12 M 36 M 



(3.16) 



7 2 



This reproduces the leading powers of Pq in the three-loop results [TT] . 



4 Heavy Quark in HQET 

Now we turn to Heavy Quark Effective Theory (HQET, see, e.g., [12111311111) and discuss 
the heavy-quark propagator. The one- loop expression for the self-energy T 1 {uj)/oj (Fig. EJ) 
in the Landau gauge with the gluon denominator raised to the power n = 1 + [L — l)e is 

, , iC F f d d k to v^v" ( k u k v 



oo 2 J ir d / 2 k-v + 00 {-k 2 ) n V ~k 2 
Using the one-loop HQET integrals 

/(^(r^r=^H^/(»,»>>. 

r/ v r(-d + m + 2n 2 )r(d/2 - n 2 ) 

nni,n 2 ) = , . , r , 

r(ni)r(n 2 ) 
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we can easily find the function F(e,u) p. 12)1 . Such functions for all off-shell HQET quan- 
tities have the same T-function structure resulting from (|4.1j) with 112 = 1 + u — e: 



F ^-( J t) ' ^ { -^%- n) Dl<n^ (4.2) 

The first T-function in the numerator, with the positive sign in front of u, comes from the 
first T-function in the numerator of 1)4.1)1 . with the negative sign in front of d, and its poles 
are UV divergences. The second T-function in the numerator, with the negative sign in 
front of u, comes from the second T-function in the numerator of 1)4.1)1 . with the positive 
sign in front of d, and its poles are IR divergences. For H(u)/u, we obtain 

N(e,u) = -2C F (3-2e). (4.3) 

If a 7^ 0, the one- loop term proportional to a should be added to E(a>). This formula 
reproduces the largest powers of (3q in the known results ^Hj at L = u/e = 1 and L = 2. 




Figure 3: One-loop heavy quark self-energy in HQET 



The heavy-quark propagator S(u) = [to — S(u;)l with the l/f3o accuracy in the 
Landau gauge is equal to 1/uj times (11.11)1 . where F(e,u) is given by (j4.2j) . ()4.3j) . Terms 
with negative powers of e in its expression via renormalized quantities form Zq. The 
anomalous dimension is given by (|1.14)1 : 



7 



N(-/3,0) 



Qf3 B(2 + f3,2 + (3)T(2 + f3)T(l-f3) 



+ 



1% 



In the general covariant gauge, the one-loop term proportional to a should be added: 



1Q 



C F — 

C F 



47T 

a* 



2a - 



1 + 1/3 



Air 



[2( 



B(2 + P,2 + P)T(2 + P)r(l-P) 
3)-8P + fp 2 + ---}. 



(4.4) 



This perturbative series converges at \P\ < 5/2. It reproduces the leading-/?o terms in the 
two- and three-loop results [T3J UHl E] . 

The renormalized expression for uS(u) is given by 1)1.15)1 . If we factor out its fi- 
dependence as in p. 16)) . then the corresponding renormalization-group invariant is given 
by fT - ^ with [IH! 



S(u) 



T(2 + u) 



2u 



-6C f 



r(-i + 2u)r(i - u) 1 

Y(2 + u) + 2^ 



(4.5) 
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(here —2a; plays the role of m). The first r-function, with the positive sign in front of 
u, produces UV renormalons, while the second one, with the negative sign, produces IR 
renormalons (Fig. EH). We can understand this from the power counting (Sect. EJ. The 
heavy-quark self-energy has a linear UV divergence which is not nullified by the Lorentz 
invariance: = 1- This is the same divergence as that of the Coulomb energy of a point 
charge in classical electrodynamics. Therefore, UV renormalons are situated at u < 1/2. 
The index of the IR divergence of the self-energy, like that of any off-shell quantity, is 
z/jr = — 2, and IR renormalons are at u > 1. 







-4 ' -3 ' -2 ' -1 ' 


[— I 2 S 1 



a 







-4 -3 -2 -1 




1 


2 " 3 " 4 



Figure 4: Renormalons in the off-shell HQET self-energy (a) and on-shell heavy-quark 
self-energy (b) 

Here we encounter a radically new situation: an UV renormalon at u > 0. It leads to 
the ambiguity AE(w)/w = (r//? )e 5//5 Aj^g/(— 2u) where r = 4CV is the residue of S(u) at 
u — 1/2. If we change the prescription for handling the pole at u — 1/2, we have to change 
the zero-energy level of HQET. Therefore, the HQET meson energy has an ambiguity 
AA = AE(w) of order A m /Po HE] (see also P2j) 

AA = -2CVe 5 / 6 ^p. (4.6) 
Po 

The structure of the leading UV renormalon at u = 1/2 can be investigated beyond 
the large-/?o limit [2D]- This approach is based on the renormalization group [HJI22]- 
renormalization-group invariant corresponding to ouS(u) is now written as 



oo 

47T 



1 + — / duS{u) exp 
Po 



o 



(3 a s (fi 



-u 



(4.7) 



instead of ljl.20j) . where the exact a s is used in the exponent, /io = —2uoe 5//6 , and 0(1/ (3%) 
is absent. The singularity of S(u) at u = 1/2 becomes a branching point 

T 

(2 

with the cut from 1/2 to +00, instead of the simple pole. The renormalon ambiguity of 
il(uj)/uj is defined, as before, as the difference of the integrals ()4.7|) below and above the 
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real axis divided by 2iri: 

. T rw 1 [ du 
AA 



fa 2m J (I_ M )^ eXP 
-2a;) exp 



47T 



2/3 r(l + a) 



2tt 



A)a s (/4))_ 



/^o« s (/io) 

47T 



(4.8) 



(Fig. El we have used T(— a)T(l + a) = — 7r/ sin7ra). But this must be just some number 
times Ajyjg, and cannot depend on uj\ 



C 



Figure 5: Integration contour 



We have to use a formula for a s (/i) more precise than the one- loop one Ijl.lOjl . The 
renormalization-group equation is solved by separation of variables: 



27r f da s 

To 



o; 



2tt 



Poa,(jJL) 2(31 
and hence 

A ms = J^exp 



47T 



2tt 



a s (n) 



log 



ft 



dlog/x. 



A 



MS 



[1 + 0(a s 



(4.9) 



(4.10) 



The UV renormalon ambiguity AA must be equal to A^jg times some number: 

AA = iVoA , A = -2C F e 5 ^^ . 

Po 

The normalization factor N is only known in the large-/?o limit: 
N = 1 + 0(l//? o ) ; 

in general, it is just some unknown number of order 1. Comparing (|4.8J) with (|4.10)l . we 
conclude that at u — > 1/2 



S{u) 



(4.11) 
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where 



0i 
ft* 



(4.12) 



The result for the power is exact; the normalization cannot be found within this approach. 

The self-energy with a kinetic-energy insertion (Fig. |HJ) can be also easily calculated 
in the large-/?o limit. In the Landau gauge, raising the gluon denominator to the power 
n= l + (L-l)e: 



(27T) 



£ & (u/) = «Cf0o 

we obtain ()4.2j) with 

iV(e,u) = 2(7 F (3-2e)V 
and hence 

AS fe (cj) = -3uAA. 



2„,i/i r 



[-k 2 ) n (k ■ v + u) 2 



g^u + 



k^k v 
^k 2 



(4.13) 



(4.14) 



This leads to the UV renormalon ambiguity of the heavy-quark field renormalization con- 
stant 



AZr 



3 AA 

2 m 



(4.15) 





Figure 6: One- loop diagrams for 

Let's now discuss the heavy-light quark current in HQET. If the light quark is massless, 
we may take | Tr of 7-matrices on the light-quark line of any diagram for T(u, 0). All 
diagrams with insertions to the gluon propagator of the one-loop diagram (Fig. Ujp), as 
well as this one-loop diagram itself in the Landau gauge, vanish due to transversality of 
the gluon propagator. Therefore, to the first order in l//3 , T(u, 0) = 1, and jj = |(7q+7 9 ) 
in the Landau gauge. This anomalous dimension is gauge- invariant, and [2] 

~ _ r a. 1 + n ( }_ 

lj ~ F 4irB(2 + !3,2 + !3)T(3 + l3)T(l-0) + \ B 2 



4tt B(2 + (3,2 + 3)T(3 + (3)T(1 - (3) \3 2 J (416) 

47T 
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UJ P 

a 



Figure 7: Proper vertex T(u,p) of a heavy-light HQET current 

from (|3.4J) and ()4.4|) . This perturbative series converges at /3o|a s | < 47r. It reproduces the 
leading- (3q terms in the two- and three-loop results [El EH HZ| • Note that jj = ^7jo (IH.9J1 
at the first order in 1//3q. 

Finally, we discuss the heavy-heavy current 

J = Zj 1 (cosh'&) J , Jo = Qjo^o, cosh i? = u • i/ 

in HQET. At one loop (Fig. |HJ), we use the Fourier transform of the Landau-gauge gluon 
propagator with the denominator raised to the power n, 

i T(d/2 — n) {2n — l)x 2 g^ u + (d— 2n)x^x u 

2 2n-l 7r d/2 r(n+l) (-X 2 + i0)d/2-n+l ' 

to obtain the coordinate-space vertex 

f l( M';coshtf) = -^^(W 

(2n - l)x 2 cosh + (d - 2n) (t + f cosh #) (f + t cosh tf) 

X 



x 2 = t 2 + t' 2 + 2tt' cosh$ . 




Figure 8: One-loop heavy-heavy vertex 
The momentum-space vertex function is expressed via the coordinate-space one as 



r(w,u/; coshtf) = y dtdt' e tujt+w,t 'T(t,t';cosh<d) . 
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Ultraviolet divergences of T(u, uj'; cosh-#) do not depend on the residual energies cu, u', 
and we may nullify them. An infrared cutoff is then necessary to avoid IR 1/e terms. 
Proceeding to the variables 



t = T 



1 + e 



t' = r 



2 ' 2 
we obtain the coefficient aAn) 



aAn) 



£j <2d—2n—2 



T{d/2 - n) 
T(n + 1) 



dr 



-d-2n-\ 



(2n - 1) coshtf(c 2 - s 2 «£ 2 ) + (d- 2n)(c 4 - s^ 2 ) 

(_ c 2 _|_ s 2^2^d/2-n+l 



where c = cosh | , s = sinh | , and the upper limit T provides an infrared cutoff. Changing 
the integration variable £ = tanh-?/;/ tanh(-$/2), we obtain 



ai(n) 



r(d/2 -n- 1) 

r(n + 1) 



-Tcosh — 
2 2 



2n+2-ci 



cosh 



2n+2-(i 



h n — 1 ) coth $ + 



d/2-n 
sinh $ 



-0/2 



cosh 2ip 



(it becomes real in the Euclidean space T 

r(i — u 



%T E ). Therefore (fTT^I . 

2u 



F(e, u) = - C F 



r(2 + w-£ 

+#/2 

dip 



cosh 2 " ^ 



-t?/2 



e 7 £jD ( £ )«A-i ( l^Tcosh- 



l—u 

(2+u- 2e) coth$ H — sinh 2-0 

smh v 



The anomalous dimension corresponding to Z-p is fll,14|) 

2(1 +/3)$cothtf + 1 



Tr 



: U F - 



3 "47r5(2 + /3,2 + /3)r(2 + /5)r(l-/5) 
In order to obtain 7 j = 7r + 7q, we add (|4.4jl : 



7J = 77^. 



1? cothf? 



3 '471-5(2 + /?, 2 + /3)r(l + /3)r(l-/3) 



4C, 



Air 



(tfcothtf- 1) . 



(4.17) 



It vanishes at 1? = as expected. It reproduces the leading (3q terms in the two-loop 
result 
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5 On-shell Heavy Quark in QCD 

Now, we turn to the on-shell mass and wave-function renormalization of a heavy quark in 
QCD at the order l//?o- It is convenient [TH] to introduce the function 

T(t) = — Tr(j> + l)E(mu(l + t)) , 
Am 

then the renormalization constants are 

Z- = 1-T(0), Z% = [1 - T'(O)]- 1 . 

At one loop (Fig. EJ), in the Landau gauge, with the gluon denominator raised to the power 
n = 1 + (L — l)e, we have 



ai(n) = —iCp 
= -2iC F 



d d k Tr(^ + 1)7^ + ft + m)Y 



fiV + 



~D7 



1 



7T 



(d-2)D 2/ \ 9N 



where 

p = mv(l + t), Di(t) = m 2 - + p) 2 , D 2 = -k 2 , 
Expanding 1/-Di (i) to the linear term in t and using the one- loop on-shell integrals 

d d k 



d/2 d-2(n 1 +n 2 



m 1 m 



m 2 - (k + mv) 2 ] ni {-k 2 ) n * 

T(d - m - 2n 2 )T(-d/2 + m + n 2 ) 



n 1: n 2 ) 



(5.1) 



M(n 1 ,n 2 



we find F(e,u) (ITT21 . 



r(ni)r(d — ni — n 2 ) 



k + mv 



Figure 9: One- loop on-shell heavy-quark self-energy 

The functions F for all on-shell quantities have a common T-function structure resulting 
from (j5.1j) with n 2 = 1 + u — e: 



\mJ 1 (o — u — e) 



(5.2) 
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The first T-f unction in the numerator, with the positive sign in front of u, comes from the 
second T-function in the numerator of 1)5.1)1 . with the negative sign in front of d, and its 
poles are UV divergences. The second T-function in the numerator, with the negative sign 
in front of u, comes from the first T-function in the numerator of 1)5.1)1 . with the positive 
sign in front of d, and its poles are IR divergences. For T(t), we obtain [2j 



N(e, u) = 2C F (3 - 2e)(l - u) [l - (1 + u - e)t] + 0(t 2 



(5.3) 



The on-shell mass renormalization constant = mo/m = 1 — T(0) with the 1 //5 
accuracy is given by ()1.11|) . ()5.2)) with N(s,u) equal to minus f)5.3|) at t — 0. Retaining 
only terms with negative powers of e, we obtain the MS mass renormalization constant Z m 
(because contains no IR divergences). Using 1)1.14)1 . we reproduce the mass anomalous 
dimension (J3.10)) . Retaining terms with e°, we get Z™/Z m (p) = m{fi)/m in the form 1)1.15)1 . 
As usual, it is convenient to express m(/z) via the renormalization-group invariant in ()1.16j) . 
Then the ratio 



(5.4) 



oo 

m (3oJ \I3qJ 



S{u) = 6C f 



r(«)r(i - 2u) 
r(3 - u) 



u 



i 

2u 



The first T-function, with the positive sign in front of u, produces UV renormalons, while 
the second one, with the negative sign, produces IR renormalons (Fig. |3Jd). We can under- 
stand this from the power counting (Sect.EJ). The QCD quark self-energy has a logarithmic 
UV divergence (iajv — 0), and hence UV renormalons are situated at u < 0. The index of 
the IR divergence of the on-shell quark self-energy is z/jr = — 1, and IR renormalons are at 
u > 1/2. 

The ratio ()5.4jl can be represented [Sj in the form ()2.4jl . At r > 1, the distribution 
function is given by the sum (J2.9j) over the UV renormalons at u = —n, n = 1, 2, 3, . . . : 



w {t) = 6C F J2 



(n + 1) (2n) 



' n\(n + 2)\ \ r 



-C F T 2 



T 



1 + 



T 



At t < 1, it is the sum ()2.8)) over the IR renormalons at u 
at u = 2: 



n + i (n = 0, 1, 2, . . . ) and 



w{t) 



C F 
C F 



T 



(2n - 1) (2n - 1)!! (2n - 5)!! r n+ ^ 



n=0 

;2-r)v^(4" 



(2n)\ 



;_4)r 



r + T 



(both of these series are easily summed using the Newton binomial expansion). Finally, 
we obtain (Fig. [10)1 



w{r) = -C F (2-t) v / t(A + t) + t 2 -66(r- V 

£ L 



(5.5) 
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At r — > 0, w(t) ~ y/r, and at r — > oo, it?(r) ~ 1/r, according to the positions of the 
nearest IR and UV renormalons. 








2 — S 




Figure 10: Virtuality distribution function 



The IR renormalon ambiguity of the on-shell mass is [TH], from the residue of S(u) at 
the leading IR renormalon u = 1/2, 



The meson mass is a measurable quantity, and must be unambiguous. In HQET, it is an 
expansion in 1/m. Its leading term, m, is a short- distance quantity - a parameter of QCD. 
The first correction, A, is a long-distance quantity, determined by the meson structure at 
the confinement scale. However, MS regularization scheme contains no strict momentum 
cutoffs. As a result, the on-shell mass m also contains a contribution from large distances, 
where perturbation theory is ill-defined. This produces the IR renormalon ambiguity (|5.fij) . 
which is suppressed by 1/m as compared to the leading term. Likewise, A contains a 
contribution from small distances, which leads to the UV renormalon ambiguity ()4.6|) . 
They compensate each other in the physical quantity - the meson mass. In other words, 
in MS the separation of the short- and long-distance contributions is ambiguous, though 
the full result is not. 

This cancellation should hold beyond the large- f3$ limit. Therefore |271] . 




(5.6) 



S{u) 



2C F N' 




«)] > 



(5.7) 




where the power is exact. The coefficients in the perturbative series 




at L 3> 1 are, according to (jl.24j) 



c n+1 = 2 1+a 2C F N^(2p ) n (l + a) n [1 + 0(1 /n)\ 



01 



a = 



2(3 { 



22 



From the Stirling formula, T(n + 1 + a) = n a n\[l + 0(l/n)], and we arrive at 



ft 



-71+ 1 



4C F AT n! (2/5 ) n (2/3 n) 2 ^o [1 + 0(l/n)] . (5.i 



This result is model-independent. 

Our calculation of T(t) also yields Zq = [1 — T'(O)] -1 at the first order in 1//%. It has 
the form (fTTH). (l5~2l with 

Nz(e, u) = -2C F (3 - 2e){l - u)(l + u - e) (5.9) 

(see (|5.3|l ). If we retain only negative powers of e, we obtain Z q {p)j Zq{^) (because Zq = 
1). Therefore, calculating the corresponding anomalous dimension by (jl.l4j) . we obtain 

lq 1Q F 4ir B(2 + (3,2 + /3)T(3 + /3)r(l - (3) \(3 2 

This difference is gauge-invariant at the level; it agrees with (|3.4j) . (|4.4jl . If we 

retain terms with e°, we get the finite combination ZqZq(^l) /Z q (fi) of the form p.!5jl : the 
corresponding renormalization-group invariant (jl. 20)1 has j2F 



= -6C F 



T(u)T(l-2u) u 2 1 
1 — it 



r(3-u) v ' 2m_ 

6 Chromomagnetic Interaction 



Now we shall discuss (22] the chromomagnetic interaction coefficient C m (/i) in the HQET 
Lagrangian. It is defined by matching the on-shell scattering amplitudes in an external 
chromomagnetic fields in QCD and HQET at the linear order in the momentum transfer 
q. All loop diagrams in HQET contain no scale and hence vanish. The QCD amplitude at 
the first order in 1/(3q is given by the L-loop diagrams with L — 1 quark loops (Fig. ITTj) . 
The results have the form (|5.2jl . The diagram of Fig.lTTk is calculated in the standard way, 
and gives 

N a (e, u) = (2C F - CU)(3 + 2u - u 2 - he + 3eu - 2eu 2 + 2e 2 - 2e 2 u) . 

We should sum the diagrams inlllb in I from to L — 1. All terms in the sum are equal, 
so that the summation just gives the factor L = u/e: 

u 

N b (e,u) =C A {h-2u-3e)-. 

In order to calculate the diagrams in Fig. fTTh. we need the triangle quark loop with the 
linear accuracy in q; it is a combination of one- loop propagator integrals (|3.1jl . Again, all 
terms in the sum in / from to L — 2 are equal, and the summation just gives the factor 
L-l: 

%T , \ _ 10 - 4m - 28e + 9eu + 23e 3 - 4e 2 u - 6e 3 (u , 

NJe.U) = —Ca ; ; 1 

v ' ; 2(1 -e) \e 
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Also, we should include the one-loop on-shell quark wave-function renormalization contri- 
bution ()5.9j) multiplied by the Born scattering amplitude (which is just 1). Finally, we 
arrive at (22] 



N(s, u) = C F N F (e, u) + C A N A (s, u) , 
N F (e, u) = 4u(l + u- 2eu) , 
2 — u — £ 

N A (e, u) = — -(2 + 3m - he - <oeu + 2e 2 + Ae 2 u) 

Z I 1 £ ) 



(6.1) 



The sum is regular at the origin e = u = 0, unlike separate contributions. It reproduces 
the leading (3 Q terms in the two-loop result 



L-l-l 





1 



L-2-l 




Figure 11: Quark scattering in an external gluon field 

Now we can easily find the anomalous dimension j m and C m (/x) with the 1//3q accuracy. 
The anomalous dimension (|1.14j) is 23 



In 



c, 



a s /3(1 + 2/?)r(5 + 2p) 

27r24(l + /?)r 3 (2 + /?)r(l -(3) 
a. 



+ o 



ft 



(6.2) 



Ca^ [l + f/?-|/3 2 + ...] . 



It reproduces the leading-/3 term of the two-loop result j23 I2H] 



7n 



A 

2tt 



(13ft - 25C A ) 



a s 
24^ 



+ 



The perturbative series ()6.2|) converges at /3 |a s | < 47r. 

The renormalizat ion-group invariant C m corresponding to C m (fi) (see (11.16)) ) has the 
form (dH]) with [23] 



S(u) 



r(«)r(i-2u) 
r(3 - u) 



4u(l + u)C F + -(2 - u)(2 + 3m)C a 



Ca 
u 



(6.3) 



The renormalon poles coincide with those in Fig. ^p. Taking the residue at the leading IR 
pole it = 1/2 and comparing with ()4.6|) . we obtain 



AO, 



7Ca\ AA 
8~C F 



m 



(6.4) 
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In physical quantities, such as the mass splitting tub* —ttlb, this IR renormalon ambiguity 
is compensated by UV renormalon ambiguities of the matrix elements in the 1/m correc- 
tion. Detailed investigation of this cancellation allows one to find the exact nature of the 
singularity of S(u) at u = 1/2: it is a branching point, a sum of three terms with different 
fractional powers of ~ — u, where the powers are known exactly, but the normalizations - 
only in the large-/9 limit. The large- L asymptotics of the perturbative series for C m can be 
found. These results have been obtained in [SB]- We shall not discuss them here, because 
they require the use of 1/m 2 terms in the HQET Lagrangian. A similar analysis of bilinear 
heavy-light currents will be presented in the next Section. 
We can rewrite C m in the form ()2.4j) with 



w(t) = C f w f (t) + C a w a (t) , 
2 + 4r + r 2 



W F {T 



w a (t 



2t 

T 

4 



vM4 + r) 
14 + 5r 



- 2 - r 



(6.5) 



9{r-l) 



Vr(4 + r) 

(Fig. these formulae can be derived in the same way as (|5.5|1 ). 









L— — 2 — 5 







Figure 12: Distribution functions w-p (dashed line) and wa (solid line) 



7 Heavy— Light Currents 

Hadronic matrix elements of QCD operators, such as quark currents j, are expanded in 
1/m 

<J> = c< ^ + i5> <6 - > + (i)' ("> 

to separate short-distance contributions - the matching coefficients C, Bi, . . . , and long- 
distance ones - HQET matrix elements <J>, <Oi>, . . . The QCD matrix element <j> 
contains no renormalon ambiguities, because the operator j has the lowest dimensionality 
in its channel. In schemes without strict separation of large and small momenta, such 



25 



as MS, this procedure artificially introduces infrared renormalon ambiguities in match- 
ing coefficients and ultraviolet renormalon ambiguities in HQET matrix elements. When 
calculating matching coefficients C, . . . , we integrate over all loop momenta, including 
small ones. Therefore, they contain, in addition to the main short- distance contributions, 
also contributions from large distances, where the perturbation theory is ill-defined. They 
produce infrared renormalon singularities, which lead to ambiguities ~ (Ay^/m) n in the 
matching coefficients C, . . . Similarly, HQET matrix elements of higher-dimensional op- 
erators <Oi>, . . . contain, in addition to the main large-distance contributions, also 
contributions from short distances, which produce several UV renormalon singularities at 
positive u. They lead to ambiguities of the order times lower-dimensional matrix ele- 
ments (e.g., <]>)■ These two kinds of renormalon ambiguities have to cancel in physical 
full QCD matrix elements <j> (JHJ) [2E| (see also 

Let's consider the leading QCD/HQET matching coefficient Cr(//) for the heavy-light 
current with a Dirac matrix T having the properties 

fT = aYf , a = ±1 , 

and (jH.fjjl . The QCD vertex function T(rnv, 0) at one loop (Fig. H3J), in the Landau gauge, 
with the gluon denominator raised to a power n, can be calculated, once and for all Dirac 
matrices, via h(d) (|3.7|) 0: 

iC F f d d k 2(d-l) + {dD 2 /m 2 + A)h-2(D 2 /m 2 + A)h 2 



2(d-l) J 7r rf / 2 D x Dl 
At the first order in l//3 , T(mv, 0) has the form ffL~TT|). with 

N(e,u) = -C F (2-u-e + 2uh-2h 2 ) . (7.2) 

In order to obtain the renormalized matrix element (Zq) X ^ T(mv , 0), we should add 
\N z (e,u)^. 

i i 
i i 

■> ' ■> — tr^ ' ■> 

7> 



V V 




Figure 13: Proper vertex T(p,p') of a heavy-light QCD current 
With the considered accuracy, all loop corrections in HQET vanish. Retaining negative 

1 /2 — 

powers of e in [ZqJ T(mv,0), we obtain Zj(fi)/Zj(fjL). The corresponding anomalous 
dimension ()1.14|) is 

„ a s 2 + (3-2(n-2-f3) 2 + (3 + 2P)(l + f3) , n f l_\ (7 ^ 

lin 7i-^ 12w S ( 2 + / 3,2 + /3)r(3 + /3)r(l-/3) + UV ' 
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in agreement with ()3.9j) and (|4.16p . Retaining e° terms, we obtain Cr(//) in the form (|1.15|) . 
The corresponding renormalization-group invariant f)l. 16|) has the form ()1.20[) with 2j 



S(u) = -C F < 



r(u)r(l - 2u) 



r(3-«) 

5 -2{n- 2) 2 
2u 



[5 - u - 3u 2 + 2ur](n - 2) - 2(n - 2) 2 ] 



(7.4) 



Comparing the residue at the leading IR renormalon u = | with AA (|4.6|) . we obtain the 
ambiguity of the matching coefficient [2] 



AC r (/i) = 3 



15 

— + r/(n-2) -2(n-2) 2 



AA 



(7.5) 



Matching coefficients for the currents with 7^° and 7^ v have identical S(u) and C* r ; 
they only differ by K 7 (a s ) in f)1.16j) . For notational convenience, we shall use the v rest 
frame. From ()7.2j) and [2j 



Ci 

C 7 HV 



C 7 o 



C* 7 HV 7 C* 7 HV 7 i 



we trivially reproduce (|3.1fij) . Taking into account [2j 

m Ci(fji) 
m(/z) C 7 o(/i) ' 

the result ()7.2j) also reproduces the corresponding formula for m/m(/i), namely ()5.3|) with 
t = 0. 

The ratio 



Ib C 7 ° 
is given by (001) with [2H] 

r(i + «)r(i - 2u) 



5(u) = 4C F - 



r(3 - tz) 



(7.6) 



It can be rewritten in the form (|2.4)l . Summing the series (|2.8jl . ()2.9)1 over the residues of 
S(u), we obtain [Hj 



tw(r) = --C F [(1 + r) vM4 + r) - r (3 + r) 



(7.7) 



(Fig. EI). 
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Figure 14: Virtuality distribution function for / Jb 



Now we shall discuss the expansion (|7.1j) for the matrix elements from B to vacuum 
of the currents with Y = 7^ c , 75 lC 7° in more detail (SOI ETJ The leading HQET current 
J— qj^ c Q has the matrix element 

<0\j(im)\B> = -i^F(fx) . (7.8) 

There are 4 dimension-4 HQET operators: 2 local ones and 2 bilocal ones. The local 
operators are full derivatives, and their matrix elements are expressed via A times ()7.8|) . 
The bilocal operators are 



O jk = i J dx T |j(0),O fc (x)} , O jm = i J dx T {j(0),O m (z)} 



(7.9) 



where O m are the kinetic operator and the chromomagnetic one in the HQET la- 
grangian. Their matrix elements are 



<Q\O jk (fi)\B> = -iy/^F(ji)G k ((j,) , <0\O jm (ji)\B> = -iy/^F(n)G m (ji) . (7.10) 
In the leading logarithmic approximation (LLA), 

<WJ /"» / , . \ A / 2/3 ° 



A 7$ 



G fc (//)=G fe -A-^-lo, 



47T 



G m (/x) — G Ti 



An 



+ 3-A 

7m0 



(7.11) 



where Gfc jm do not depend on y^, and y°' m are mixing anomalous dimensions of Ojkj m with 
local operators, see [HI]. The renormalization-group invariant QCD matrix elements of the 
pseudoscalar current and the axial one are in LLA 



+ 



a s {m) 
An 

1 

2m 



MA 



2 A, 



47T 



7m0 



A + G fe + G r . 



An 



'mil 

2/3 



f7.12) 
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The next-to-leading corrections are also known [3~2~1 131 j . 

We are interested in UV renormalon ambiguities of the matrix elements of Ojkjm- By 
dimensional analysis, they are proportional to AA times the matrix element of the lower- 
dimensional operator J with the same external states. We may use quark states instead 
of hadron ones. Specifically, we consider transition from an off-shell heavy quark with 
residual energy u < to a light quark with zero momentum, this is enough to ensures the 
absence of IR divergences. For Ojj., all loop corrections to the vertex function vanish. The 
kinetic-energy vertices contain no Dirac matrices, and we may take \ Tr on the light-quark 
line; this yields k a at the vertex, and the gluon propagator with insertions is transverse. 
There is one more contribution: the matrix element of J should be multiplied by the heavy- 
quark wave-function renormalization Zq 2 , which has an UV renormalon ambiguity (|4.15|) . 
Therefore 



AG k (^) = -^AA (7.13) 

(this ambiguity is /i-independent at the first order in 1/ (3q). 

k 





uo kTTuv 1 ~k~ > 



Figure 15: Matrix element of Ojk 

For Oj m , a straightforward calculation of the diagram in Fig. [TBI gives the bare matrix 
element of the usual form (|1.11|) . (j4.2j) with 

N(e,u) = -6C F C° m -. (7.14) 
m 

The renormalization-group invariant matrix element has the form (jl.20|) with fi = —2cue~ 5 ^ 6 
and 

c( \ rr r f , u f T(-l + 2u)T(l - u) 1\ 

S(u) = -6C F CM~ { f ^ + Yu ) • (7.15) 

Taking the residue at the pole u = 1/2, we find the UV renormalon ambiguity C m (/io) AA/m 
times the matrix element of J, and we obtain |2Ej 

AG m (fi) = 2AA (7.16) 

(again, //-independent at this order). 

In the full QCD matrix elements (J7.12j) . the IR renormalon ambiguities ()7.5j) of the 
leading matching coefficients Cp are compensated, at the l//?o order, by the UV renormalon 
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Figure 16: Matrix element of 0j r 



ambiguities of the subleading matrix elements AA ()4.6|) and AGk, m (|7.13jl . (|7.16jl . This 
cancellation must hold beyond the large-/?o limit. The subleading matrix elements are 
controlled by the renormalization group. The requirement of cancellation allows one to 
investigate the structure of the leading IR renormalon singularity of Cr |31j . 
In the large-/3 limit, 



AG* 



■*AA 
2 



AG r 



ImO 



AA 



see f!7.13j) . ()7.1(ij) . (|7.11jl . In general, they must be equal to times some numbers: 



AGi 



■-iVxAo, AG m = N 2 2 



To 



A r 



(7.17) 



2 V 7m0 , 

(see (I4.1()jl ). The normalization factors Ni^ are only known in the large-/3o limit: 
N t = 1 + 0(1/ (3 ) ; 

in general, they are just some unknown numbers of order 1. Using ([4.9)1 . we can represent 
the UV renormalon ambiguities of the 1/m corrections in ()7.12|) as exp[— 27r/(/5 a s (/io))] 
times a sum of terms with different fractional powers of a s (//o)/(47r). It is convenient to 
replace log[a s (/io)/(47r)] — > [(a s (/io)/(47r))' 5 — l]/6, and take the limit 5 — > at the end of 
calculation. 

In order to cancel this ambiguity, we should have the branching point 



Sr(u) 



instead of a simple pole ()7.4|) . Then (see (|4.8|) 1 



AC r = ~r exp 

Po 



2tt 



/3 a s (/i ) 



+ ^ 



47T 



(7.19) 



The requirement of cancellation of the ambiguities in ()7.12|1 gives for T = 1, 7 



S r (u) 



(i-Ti) ^ 



7o_ 

2/3 



iogi^-^i+ A 







±1 + 



7m0 



7m0 



^2 (| - «) ^ 



(7.20) 
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The similarly requirement gives for T = 7*, 7*7° 



Sru) 



1+ 



0i 
-■vi 



To 



7m0 



±1 



To! 

7m0 



(7.2i; 



Here 



jv[ = mr ( 1 + A ) flj* , jv* = jv a r [ 1 + 



A 7m0 



2$ 2/5 



01 TVn.0 
2^ 2/3 



sec 

flUOD). Corrections (§ - u) were calculated in [31]. In the large-/?o limit, the simple 
pole behaviour with ()7.5|) is reproduced. 

The asymptotics of the perturbative coefficients c L at L ^> 1 is determined by the 
renormalon singularity closest to the origin. Similarly to (|5.8j) . we obtain, for T = 1, 7 , 



0i 



2C F n\ (2(3 ) n (2(3 n) 2 ^ 



To 



76" 



-^-log2/? n±l + 

2po 7m0 



iVn 



-JVi+ 2-J2- ) JV 2 (2A)n)~ a * 



7m0 



(7.22) 



and for T = 7*, 7*7° 



ft 



-'ra+l 



2C F n! (2(3 ) n (2(3 n) 2/3 S 



7o_ 
2/?o 



log 2(3qu + ^ I ±1 



7m0 



2 3 



7m0 



, TWO 

N 2 (2/3on) 2 "o 



(7.23) 



Corrections 0(l/n) were calculated in [3*T] . 

For the ratio /#* the Borel image of the perturbative series is 



S(u) 



3 n v 
(2 ) 



0i 



7m0 



7m0 



, , 1 . w 

^2 (| - ^ 



(7.24) 



and the asymptotics of the coefficients is 



c n+1 = j|c F n! (2/3 ) n (2/3 n)% 



7m0 



7m0 



N 2 (2fi n) 2 ™o 



(7.25) 



31 



8 Heavy— Heavy Currents 

First we consider the leading matching coefficients for the currents cTb at $ = at 1//3q 
order [HB]. The one-loop vertex T(mbV,m c v) (Fig. I17jl contains integrals 

d d k 



[m 2 b - (k + m b v) 2 } ni [m 2 c - (k + m c v) 2 f 2 (-k 2 ) n ' 

The denominators are linearly dependent. We can multiply the integrand by 

rrib [m 2 c — (k + m c v) 2 } — m c [ml — (k + m b v ) 2 ] 
(m b -m c )(-k 2 ) 

thus lowering n\ or ri2, until one of these denominators disappear. Remaining integrals are 
single-mass (|5.1|) . We have 

, , r(rf-2n-l)r(-rf/2 + n + l)m^~ 2n - 1 $(m c /m 6 )-mf- 2n - 1 $(m 6 /m c ) 

T(d-n-l) m b -m c , g ^ 

d-1 2/i 2 + (d-2n-2)/i-d + n + l 

$(r) = r H . 

y ' d-2n-3 d-n-1 



i 

:> * =» , 

a 




Figure 17: Proper vertex p') of a heavy-heavy QCD current 
Adding the on-shell wave-function renormalization ()5.9|) for b and c, we obtain 



F(e,u) 



Vm 6 m c y r(3-w-e) V ; V 



iV(e, u) = 2Cj 



(n - 2) 2 - wr^n - 2) + 2e:(n - 2) - uerj 



I — 2u 

+ (3 - 2e)- — (1 - u - m 2 + M£)i? c 

1 + 2-u 



w), (8.2) 
4 — u 2 — As — 2ue 2 ' 



l + 2u 



Ri 

8.3) 



where 



L 



smh ^ — 5-^— 

i?o = cosh — , R x = — — f — 

I smh ~ 



L = lot 



m b 



m r 



for the on-shell QCD matrix element (which is equal to the matching coefficient, because 
all loop corrections in HQET vanish). The corresponding anomalous dimension ()1.14|) 
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reproduces jj n (|3.9|) . because jj = at $ = 0. The function ()1.2H) for the matching 
coefficient is 



S( u ) = C F ^2 
+ 3 



r(u)r(i-2u) 



r(3-u) 

(l - 2u)(l - u - u 2 ) 
1 + 2u 



(n - 2) 2 - «ry(ri - 2) 



4-w 2 
l + 2w 



i?0 



(n- l)(n-3) 



u 



U) 



with 



/i = e 5/6 y/m b m c 

(see (11.19)) ). There is no pole at it = 1/2, the leading IR renormalon is at u = 1. Therefore, 
the IR renormalon ambiguity of the matching coefficients at i? = is ~ (A^/m^c) 2 //^. 
For the vector and axial currents 



6C f 



T(u)T(l-2u)l-u-u 2 
r(3 -u) 1 + 2u 



-Rx + (1 - 2u)R c 



S. 



u) = 2Ci 



L , L 
— coth — 
2 2 

rMrri - 2u) 



l i 



(i + 2«)r(3 - «) 

3L coth 8—1 — L coth 



-(3 - u + u 2 )Ri + (1 - 2u)(l - it - M 2 )i? c 



15) 



6 u + 



The matching coefficients don't depend on /x, //, and are given by (jl.2()j) : 



47T 

a s (^o) 



H^ = C F ^ 



coth ■ 
2 2 

3L coth 

2 
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-Lcoth^-G j 



5.6) 



Now we consider the general case 7^ [2B]. For a generic Dirac matrix V satisfying 

T „iy = 2/i(d)r, 

there are 4 leading HQET currents in the expansion: 



J = E^ + i E G.s. + ^ E g',o; + on/ 

i j i 



b,c) > 



1.7) 



The one- loop vertex r(m;,t;, m c t>') in the Landau gauge, with the gluon denominator raised 
to a power n, is 

d d k 7^ + m c ^' + m c )r(^ + + mj ) ) r y fl 



(27r) d (-k 2 ) n (-k 2 -2m b v k)(-k 2 - 2m c v> ■ k) ' 
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Using Feynman parametrization, we have 



iC F n(n + 1) f (1 - x - x'Y~ l dx dx' d d k' 
Ti = — I - — , — — x 



n d/2 J ( fl 2 _ fc /2)n+2 

7jt($' + ?n c (l — a/)^' — m h xf + m c )r(^' + m b {l — x)f — m c x'f + m h )^ , 
k' = k + m;,xf + m c x'v' , 

a 2 = (rribxv + m c x'v') 2 = mix 2 + mix' 2 + 2m b m c xx cosh?? . 

We calculate the loop integral and substitute x = £(1 + z)/2, x' = £(1 — z)/2, 

2 ,2 ,i±i? . , £,=bi? 
a = mbTn c £ a + a_ , a± = cosh — h z smh — - — . 

The £ integration is trivial, and we obtain Hi having the form (jl.llj) . without the leading 
1 for -^2,3,4; the functions Fi(e,u) have the form ()8.2|) with 



Ni(e,u) = C F < 



dz 



N 4 (e, u) 



(a + a_) 1+u 



(I _ r 2 U 

a + a_ h 2 (l - 2u) - - —u(l - 2u) 



+ (1 + u — e) (2 — u — e) cosh $ + (cosh L + z sinh L)u(2 — u — e) 

(3 - 2e)(l - + u - e)(r u + r~") 
+i 



No(e,u) = -C f 



u 



dz 



N 3 {e, u)= -C F 



2 J (a + a_y+ u 
-l 

+i 

u f dz 



e L (l + z) 2 h 



(l-2u) + (l-z)(2-u-e) 



2 J (a + a_y+ u 
-l 



e~ L {l-z) 2 h 



(1 - 2m) + (l + z)(2-u-e) 



hu(l-2u) fdz(l-z 2 ) 
(a + a_y+ u 



At i? = 0, the integrals can be easily calculated, and iVi + iV2 + ./V3 + ./V4 reproduces ()8.3j) . 
The anomalous dimension ()1.14|) corresponding to ()8.8|) is 7 jn — 77, see (J3.9|) . ()4.17jl . The 
functions lfT2TJ) 

r( M )r(i-2 M ) Ar/n iv,(o,o) 

= r( 3 - M ) - 

have the leading IR renormalon pole at u = 1/2, thus producing the ambiguities (j2.2j) 
^(0,1/2) AA 



AE = -- 



3Cp ^mb m c 
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in the matching coefficients. It is easy to calculate A^O, 1/2) using the integrals 



+1 +i , 

/dz 4 cosh -j f zdz 4sinh-| 

(a + a_) 3 / 2 = coshtf + 1 ' J (a + a_) 3 / 2 = ~coshtf + 1 ' 

-l -l 

We obtain |2E| 

Vcosht' + l AJ \m c m b J 

1 AA 1 AA 

A# 2 = — r~s , = — r- n , &H 4 = . 

cosh 2m c cosh v + 1 2m;, 

They do not depend on the Dirac matrix V in the current. As expected, AH 4 + AH 2 + AH 3 
vanishes at i? = 0. 

In matrix elements of QCD currents (|8.7jl . these IR renormalon ambiguities in the 
leading matching coefficients Hi must be compensated by UV renormalon ambiguities in 
matrix elements of the sub leading operators Oj, 0[. There are two kinds of subleading 
operators - local and bilocal. First we consider local operators, whose coefficients are 
completely fixed by reparametrization invariance 



J2 G l O i = c v , r + H 2 fY + H 3 Vf + H 4 j>Tf) %Jp + 2iip (H 2 T + H 4 Tf) 

(8.10) 



local 



+ 2 (H[ T + H' 2 i)Y + H' 3 Tf + H' 4 jTf) iv' ■ D 



where H[ are derivatives in the argument coshi/ 1 , and similarly for J^G^O-. These local 
operators contain either or D ll . Let's decompose these derivatives into components in 
the (v, v') plane and those orthogonal to this plane. Projection of _D M onto the longitudinal 
plane is 



D 



(Vp v ' ■ D + v'^ v ■ D) cosh i? — v • D — v' v ' • D 
sinh 2 1? 



and similarly for D ^. All operators with longitudinal derivatives can be rewritten, using 
equations of motion, as full derivatives of the leading currents Jj. When we are interested 
in matrix elements from a ground-state meson into a ground-state meson, we may replace 

id^Ji —> A(v - v')^Ji . 

In this case, projecting onto the longitudinal plane means 

- v u cosh i? — v' <— - v' cosh # — v u 

iD —> A— -, -iD — > A— -. 

M cosh •& + 1 ' M cosh d + 1 
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The longitudinal part of the local contribution (|8.10|) to the QCD matrix element 

<J> is easily derived by this substitution. It, clearly, has an UV renormalon ambiguity 
proportional to AA. Matrix elements of operators with transverse derivatives cannot be 
written as matrix elements of the leading currents Ji times some scalar factors, they require 
new independent form factors. Therefore, they contain no UV renormalon ambiguities, 
which should have the form AA times lower-dimensional matrix elements of Jj. The above 
derivation is exact (not only valid in the large-/?o limit). At the first order in 1//3q, we may 
replace H\ — > 1, -£^2,3,4 — ► 0, H[ — > 0. The contribution of the local sub leading operators 
to the ambiguity of <J> in this approximation is [2H] 



1 - 



1 



1 1 \ AA < J > 1 f AA <j- AA <J- 

Trie rrib J 2 1 coshi^ + l \2m c 2 



cosh r) + 1 J V m c rrib J 2 ' cosh d + 1 V 2m c * ' 2rrib 
Now we turn to bilocal subleading operators, and consider the operator 



darT{ji(0),O te (x)} 



with the insertion of the c-quark kinetic energy. It appears in the expansion ()8.7j) with 
the coefficient Hi. The one-loop vertex (Fig. fTHj) with the gluon denominator raised to the 
power n is 



ai{n) 



.C F 
2m r 



d d k f 1JX kW 



k ■ V' + LU'J (k ■ V + Lj)(k ■ V' + LU')(-k 2 y 



n [g»u+ _ k2 



where k± = k — (k ■ v') v'. We are interested in the UV renormalon at u = 1/2; therefore, 
to make subsequent formulae shorter, we shall calculate F(u) (|1.18jl instead of the full 
function F(e,u), and omit terms regular at u = 1/2. We also set u' = to for simplicity, 
and obtain 



F{u) 



— 1-!t—u{—2uj) 



2u 



+ 



2m c 
d 4 k 



d 4 k 



1 



7T 2 (-jfc ■ V' ~ Lu)(-k 2 Y+« 

cosh'i? 



7T 2 (-jfc • V - L0)(-k ■ V' - L0) 2 (-k 2 ) u 



+ 



where dots mean integrals without linear UV divergences at u = (and hence having no UV 
renormalon singularity at u — 1/2), and —2uj plays the role of m in the definition (j!.18|) . 
The first integral is trivial (see (j4.1j0 : 



--A-2UJ) 



-l+2u 



7T 



d A k 



r(-i + 2«)r(i 



u 



(-fc • v' - uj)(-k 2 y+ u 



r 1 



u 
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Figure 18: Kinetic-energy insertions into the c-quark line 

For the second one, we use the HQET Feynman parametrization [15J: 

d A k 



-l+2u 



7T Z 



(-k-v- u)(-k ■ v' - u) 2 {-k 2 ) u 

y' dy dy' d 4 k 



i 1 ( 2^r i+2 " r(3 + M) 

vr 21 T(u) J [-tf - 2yv k - 2y'v' ■ k ~ 2u(y + y')f +u 



iu{-2u) 



-l+2u 



y' dy dy' 



[y 2 + y' 2 + 2yy>cosh$ - 2u{y + y')} 1+u ' 
The substitution y = (-2u;)f (1 - z)/2, y' = (-2w)£(l + z)/2 gives 

i x - u di{l + z) dz 



I = 2u 



l+u 



[(cosh 2 f -z 2 sinh 2 f)£+l] 
Then the substitution (cosh 2 | — z 2 sinh 2 | ) £ = r\ leads to the factored form 



I = 2u 



oo +1 

rj 1 ~ u dri 



dz 



j (7j + i)i+«y [cosh 2 !-^ siri h 2 f] 



2 #1 2-« ' 



where 



r] l - u di] r(-l + 2tz)r(2 - u) 



(v + i) 



l+u 



r(i + u) 



Collecting all contributions, we obtain (jl.21j) 
-2uj r(-l + 2u)r(l - u 



5(u) = C f 



r(i + u) 



X 



1-J. 

1 + u(l — u) cosh-# J 



dz 



/ [cosh 2 f -z2 sinh 2 f] 
where dots mean terms regular at u = 1/2. 



2-u 



(8.12) 
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The residue at the pole u = 1/2 can be obtained using 



+1 

r dz 1 

I [cosh 2 f -z2 sin h 2 f] 3/2 ~ coshtf + 1 • 

The corresponding UV renormalon ambiguity is given by (J2.2j) with — 2uj instead of m. 
Adding also the external-line renormalization effect AZ c /2 (|4.15|) . we obtain 

f f \ AA 



2 cosh i? + 1 / 2m c 

The contribution of the bilocal operator with the 6-quark kinetic energy contains my, instead 
of m c . Therefore, contribution of all bilocal operators with kinetic-energy insertions to the 
ambiguity of <J> is 



1 1 \ / 1 1 \ AA ~ 

2 _ coshtf + 1 ){m c + m b ) ~ <Jl> ' (8 ' 13) 

Matrix elements of bilocal operators with a c or 6-quark chromomagnetic insertion cannot 
be represented as matrix elements of the leading currents Jj times scalar factors - they 
require new independent form factors. Therefore, they have no UV renormalon ambiguities, 
which should be equal to AA times lower- dimensional matrix elements. 

Summing the ambiguity AHi<Ji> ()8.9j) of the QCD matrix element <J> due to the 
IR renormalon in the leading matching coefficients Hi at u = 1/2, the ambiguity 1)8.11)1 due 
to the UV renormalon in the local subleading operators at u = 1/2, and the contribution 
of the bilocal subleading operators (|8.13|) . we see that they cancel, at the first order in 
1//3q, for any Dirac structure T of the current J |2l 
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